STATIONARY AND NONEQUILIBRIUM FLUCTUATIONS IN 
BOUNDARY DRIVEN EXCLUSION PROCESSES 



C. LANDIM, A. MILANES, S. OLLA 

Abstract. We prove nonequilibrium fluctuations for the boundary driven 
symmetric simple exclusion process. We deduce from this result the stationary 
fluctuations. 



1. Introduction 

In the last years there has been considerable progress in understanding stationary 
non equilibrium states (SNS): reversible systems in contact with reservoirs imposing 
a gradient on the conserved quantities of the system. In particular, large deviation 
properties has been studied for boundary driven one-dimensional symmetric simple 
exclusion processes and references therein). 

One of the most striking typical property of SNS is the presence of long-range 
correlations. For the symmetric simple exclusion this was already shown by H. 
Spohn in the pioneering paper jl(J|. But we noticed that a mathematical proof of 
the convergence of the fluctuation fields to the corresponding Gaussian field was 
missing from the literature. The purpose of this paper is to fill this gap. 

We consider the symmetric exclusion process in an open lattice of length N. 
Particles jumps to nearest neighbors performing simple symmetric random walks 
with the exclusion rule: a jump is suppressed if site is already occupied. At the left 
boundary particles are created with rate a and annihilated with rate 1 — a. On the 
right boundary this is done with rates (3 and 1 — (3. To keep notation simple we 
restrict to the one dimensional nearest neighbor case. Extensions to more dimension 
and more general jumps rates is straightforward (see remarks 12 . 21 and \'2 . 3|) . 

If a = (3 = p, the Bernoulli product measure with probability p is stationary 
and reversible for the dynamics. But when a ^ (3 the stationary measure has 
correlations and is not explicitly computable. We denote by rj(x) = or 1 the 
occupation variable of site x. It is easy to prove that < ri([Nu\) > ss —* p(u) — 
((3 — a)u + a. 

The fluctuation field is formally defined as the random distribution on [0, 1] 

1 

Y N (u) = */*0fa(*) ~ P{u)) ■ (1.1) 
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We prove in this paper that, under the stationary measure for the process, Y 
converges in law to the centered Gaussian field Y on [0, 1] with covariance 

< Y{u)Y(v) > = X {P(u))5(u -v) - {f3- a) 2 (-A)-\u, v) , (1.2) 

where A is the Laplace operator with Dirichlet boundary conditions, and \(p) — 
p(l — p). In the 1-dimensional case we have more explicitly (— A) _1 (u,u) = u(l — v). 

The strategy we use to prove this result is to study first the convergence of the 
nonstationary fluctuations. If we start with some non-equilibrium density profile 
< T]o([Nu}) >= p(Q, u), then at the diffusive time scale we have < r] N 2 t ([Nu]) >^ 
p(t,u), where p(t,u) is solution of the heat equation with initial condition p(0,u). 

We then consider the time-dependent fluctuation field 

1 N ^ 

Y N (t, u) = -= J2 S(u - x/N)( mH (x) - p(t, u)). (1.3) 

V X— 1 

The main point of the proof is to show the convergence of Y (u, t) to the solution 
of the stochastic linear partial differential equation 

d t Y(t, u) = AY(t, u) - V (y/2 X (p{t, tO) W(t, «)) , (1.4) 

where W(t,u) is the standard space-time white noise. If we start in the stationary 
state, p(t, u) = p(u) for all t. In this case the distribution valued process Y(t, u) is a 
stationary Gaussian process and its invariant distribution is given by the Gaussian 
field Y with covariance given by ljl.2|l . 

This article presents a rigorous proof of the results described above and presented 
in 0] . Article 0] also contains the connection between the large deviations and the 
small fluctuations proved here, showing that the inverse of the covariance <|1.2t is 
given by the second functional derivative of the large deviations rate function. 

2. Notation and results 

For N > 1, let A N = {1, . . . ,N - 1}. Fix < a < (3 < 1 and consider the 
boundary driven symmetric simple exclusion process associated to a, (3. This is the 
Markov process on {0, 1} Ajv whose generator Lm is given by 

AT-2 

(L N f)( V ) = "£{f(**' x+i v)-m} 

x=l 

+ [a[l - rj(l)} + (1 - a)»;(l)}{/(«7 1 T 7 ) - f(rj)} 

+ {/3[1 - r,(JV - 1)] + (1 - (3) V (N l)}^-^) - /(r?)} . 

In this formula, rj — {rj(x), x e Ajv} is a configuration of the state space {0, 1} A ™ 
so that rj{x) = if and only if site x is vacant for 77; a x,v rj is the configuration 
obtained from r\ by interchanging the occupation variables 7](x), rj{y): 

!rj(z) \iz^x,y, 
r](y) if z = x , 
T)(x) tfz = y; 

and <j x rj is the configuration obtained from rj by flipping the variable rj(x): 
, x nVv s / V(z) if z ^ x , 

(CT ^ )(z) - ( i-r,{z) ifz = x. 
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Hence, at rate a (resp. 1 — a) a particle is created (resp. removed) at the boundary 
site 1 if this site is vacant (resp. occupied). The same phenomenon occurs at the 
boundary x — N — 1 with in place of a. 

This finite state Markov process is irreducible and has therefore a unique station- 
ary measure, denoted by v£ a . For < 7 < 1 , denote by vOf the Bernoulli product 
measure on {0, 1} Ajv with density 7. If a = (3, an elementary computation shows 
that v£ is the invariant measure and that the process is reversible with respect to 
this stationary state. On the other hand, if a < j3, it is known since ^U] that the 
invariant state has long range correlations. 

Static picture. For N > 1, denote by 71-^ the measure on [0, 1] obtained by 
assigning mass iV" 1 to each particle: 

n N ( V ) = AT 1 £ r,(x)S x/N , (2.1) 

xGAjv 

where S u is the Dirac measure on a. It has been proved in |S] that under the 
stationary state v& b t ne empirical measure ir N converges to the unique solution of 
the elliptic equation 

Ap = 0, 

P (0) = a , p(l) = (3 . 
We denote the solution of this equation by p — p a ,p- 

Once a law of large number has been proved for the empirical measure under the 
stationary state, it is natural to consider the fluctuations around the limit. Let 

p N (x)=E <0 [ V (x)}. 

Since E v n [LnT)(x)] = for all 1 < x < N — 1, an elementary computation shows 
that p N is the solution of 

' (A NP N )(x) = for 1 < x < N - 1 , 
p N (0) = a, p N (N)=f3, [Z - Z) 

where A^v is the discrete Laplacian: (A N H)(x) = N 2 {H(x+l)+H(x-l)-2H(x)}. 

In the case of the symmetric simple exclusion process, p N (■) is just the linear 
interpolation between ^(0) = a, p N {N) = (3. 

To define the space in which the fluctuations take place, denote by Cq([0, 1]) the 
space of twice continuously differentiable functions on (0, 1) which are continuous 
on [0, 1] and which vanish at the boundary. Let —A be the positive operator, 
essentially self-adjoint on L 2 [0, 1], defined by 

-A -— 

dx 2 ' 

V(-A) = C 2 ([0,1]). 

Its eigenvalues and corresponding (normalized) eigenfunctions have the form X n = 
(rnr) 2 and e n {u) = \/2 sin(n7rw) respectively, for any n € N. By the Sturm-Liouville 
theory, {e„, n € N} forms an orthonormal basis of L 2 [0, 1]. 

We denote with the same symbol the closure of —A in L 2 [0, 1]. For any nonneg- 
ative integer k, we define the Hilbert spaces Hk — T>({—A} k / 2 ), with inner product 
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(f,9)k = ({-A} fc / 2 /, {-A} fc/2 .g), where (-,-) is the inner product in L 2 [0, 1]. By 
the spectral theorem for self-adjoint operators, 

+00 

H k = {/eL 2 [0,l]: 5> 2fc Cf, e „) 2 <oc} , 

n=l 

+00 

U,g)k = ^2{nn) 2k (f,e n )(g,e n ) . 

71=1 

Moreover, if H-k denotes the topological dual space of Hk, 

+00 

H- k = {/GP'(0,1): ^n- 2fe (/,e„) 2 < 00}, 

+ 00 

{f,9)-k = ^2{nn)- 2k {f,e n )(g,e n ), 

71=1 

where (/, •) represents the action of the distribution / over [0, 1] on test functions. 
Fix k > 5/2 and define the density field Y N on H-k by 

Y N {H) = TV- 1 / 2 £ H(x/N){ V (x)-p N (x)} . (2.3) 

xeAjv 

For fc > 1, denote by q the Gaussian probability measure on H-k with zero mean 
and covariance given by 



E q [Y(H)Y(G)} = (2.4) 
dux{p{u))H{u)G{u) - (p-a) 2 / du[(-A)~ 1 H(u))G(u) . 



Theorem 2.1. Fix k > 5/2 and denote by the probability measure on H-k 
induced by the density field Y N defined in l|2.3|l and the stationary measure a- 
As N I 00, qN converges to q. 

This result follows from Proposition 13. 51 which is proved in Sectional 

Remark 2.2. The same statement holds in higher dimensions for a symmetric 
exclusion on the set An x T^ 1 , where T^r is the discrete torus of length N . The 
dynamics is periodic in the d—1 directions orthogonal to the gradient of the density. 
In this case the space correlations are given by Q2.4fl . where A is the Laplacian 
with periodic boundary conditions in the (d-1) dimensions and Dirichlet boundary 
conditions in the first coordinate. The proof is an elementary extension of the 
one- dimensional case. 

Remark 2.3. We may also consider a boundary driven symmetric simple exclusion 
process in which the occupation variables rj(x), n(x + y) are exchange at rate p(y) 
for a finite range irreducible probability p{-) . In this case, the Laplacian is replaced 

by the operator Y,i,j=i "■.,<>„ <>„ ■ where a it j = Y^yV^y^Piv)- 
3. NONEQUILIBRIUM FLUCTUATIONS 

We prove in this section the dynamical nonequilibrium fluctuations of the bound- 
ary driven exclusion process. We start with the law of large numbers. 
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Fix a density profile po : [0,1] — > [0,1]. Consider a sequence {/x , N > 1} 
of probability measures on {0, l} Ajv such that for every continuous test function 
H : [0, 1] -> R and every e > 0, 



Jim p N {n N {H)- J H(u) Po (u)du >e} 



= 



Denote by P^n the probability on the path space _D(R + ,{0, l} Ajv ) induced by 
the Markov process with generator Ln and the initial measure p N . Denote by 
tt^ the empirical measure associated to the state of the process at time t: tt^ = 
N~ x T]N2 t (x)5 x /N- It follows from the usual hydrodynamic limits techniques, 
adapted to the boundary driven context (cf. sections 4 and 5 in |7], and [S]) that 
for every t > 0, every continuous test function H : [0, 1] — > R and every e > 0, 



| 7rf (H) - j H{u)p{t,u)du > e} 



lim 1 

where p(t, u) is the unique solution of the heat equation 

d t p = Ap , 

p(0, •)=/*>(•). 

p(',0) = a, p(-,l)=/3- 
Furthermore is valid the following replacement lemma: 



= 



(3-1) 



Lemma 3.1. Let ^(n) a local function, and ^(p) = E v (ty), where v p is the 
Bernoulli measure with probability p. Let G(s,u) a continuous function on M + x 



[0,1]. Then 



( f T 1 JV_1 x f 1 
lim sup E I ds —Y] G 0> T T )^(r x VN 2 s) - / 



G(s,u)^(p(s,u)) du 



= 



The proof is given in chapter 5 of pj| , adapted to the open boundary situation. 

We now turn to the fluctuations. Consider a sequence {p, N : N > 1} of proba- 
bility measures on {0, l} Ajv . Let 

p N (x) = E^n [t]{x)] , <p N (x,y) = E^N[v{x)\n(y)] , 

for x, y in An, x < y. In this formula, E^ [/;<?] stands for the covariance of / and 
g: E^[f;g] = E^fg] ~ E^f] E^g]. We extend the definition of p N and ip N to the 
boundary of Ajv by setting 

p N (Q) = a, p N (N) = f3, <p N {x,y) = 

if x or y does not belong to Ajy. Assume that there exists a finite constant Co such 
that 

sup N\p N (x + l)- p N {x)\ < C , N max \<p N (x, y)\ < C . (3.2) 

0<x<N-l x,yeA N 

x<y 

Assume furthermore that p N converges weakly to a profile po in the sense that for 
every continuous function H : [0, 1] — * M, 



E H{x/N)p N {x) = [ duH 



(u)po(u) 



(3.3) 
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It follows from assumptions 1(3. 2 jl . ((3.3(1 and from Chebyshev inequality that 
under (j, the empirical measure ir N , defined in l|2.1(l . converges to po(u)du: For 
every S > and every continuous function H : [0, 1] — > K, 

lim p N { n N (H) - [ duH{u) P() {u) > s) = . 

N^cc I Jo J 

In particular, by the law of large numbers stated in the beginning of this section, 
for every t > 0, the empirical measure ir^ converges to the absolutely continuous 
measure whose density is the solution of the heat equation ((3.1(1 . 

We prove in 1(4.2(1 that the stationary state v£ a satisfies the assumptions 1(3.2(1 . 
(13. 3(1 . It also easy to verify that this property is shared by product measures 
associated to Lipschitz profiles. 

Let pf be the solution of the semidiscrete heat equation 

d s p™{x) = A N p s (x) , x e Ajv ; 

p»(x)=p N (x), xeA N ; (3.4) 
p?(p)=a, p?(N)=(3, s>0. 

Fix k > 5/2 and denote by Y t N the density fluctuation field which acts on smooth 
functions H in Ti.k as 

Y t N (H) = N-^ H(x/N){ Vm ,(x)-p?(x)}. 

x£iA N 

Notice that time has been speeded up by N 2 . Denote by Qm the probability 
measure on D([0,T],Ti.-k) induced by the density fluctuation field Y N introduced 
above and the probability measure (j, . 

Assumptions Q3.2p. (|3.3|l ensure tightness of the sequence Qn and permit to 
describe the asymptotic evolution of the field Y as the sum of two uncorrelated 
pieces: a deterministic part characterized by the heat kernel and a martingale. 
This is the content of the first result. Denote by {T s : s > 0} the semigroup 
associated to the operator A. 

Proposition 3.2. Fix T > and a positive integer k > 5/2. The sequence Qn is 
tight on D([0,T],H-k) with respect to to the uniform topology. All limit points Q* 
are concentrated on paths Y t such that 

Y t (H) = Y Q {T t H) + W t (H) , (3.5) 

where Wt(H) is a zero-mean Gaussian variable with variance given by 

t ,i 

ds / du X (p( S ,u))(T t „ s VH) 2 (u) , 
o Jo 

and p s is the solution of the heat equation (13.1(1 . Moreover, Yq and Wt are uncor- 
related in the sense that Eq* \Y$(H)Wt(G)] = for all functions H, G in Cq ([0, 1]) 
and all < t < T . 

Proof. The proof of tightness of the sequence Qn is left to the end of this section. 
To check the properties of the limit points, fix a smooth function H in Cq([0, 1]). 
An elementary computation shows that 

{d t + N 2 L N )Y t N {H) = Y t N (A N H). 
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Observe that no boundary term appears in the right hand side of the above equation. 
In particular, defining 

rf (H) = N 2 {L N Y S N (H) 2 - 2Y S N (H)L N Y S N (H)} , (3.6) 

it follows that 

Mt> N {H) = Y t N (H) - Y N (H) - f dsY s N {& N H) , (3.7) 

Jo 

M?< N (H) = {Ml< N {H)} 2 - fdsT^H), 

Jo 

are martingales. A simple computation shows that 

N-2 

r s( H ) = ^£[^ 5 (*+i)-^*)H(v;viO(^)} 2 ( 3 - 8 ) 

x=l 

+ N-^WnH^N - 1)/N) 2 { VN 2 S (N - 1) - (3} 2 , 

where Vat H(x/N) stands for the discrete derivative: (V mH)(u) = N {^[(u+N^ 1 ) — 
H(u)}. Observe that the last two boundary term on the above equations are of order 
N- 1 . By Lemma ITT1 as N | oo, rf (H) converges to 2 J* x{p{-% u)){VH{u)) 2 du, 
where p is the solution of the heat equation (|3.1|) . 

Fix a limit point Q* of the sequence Qn- It follows from l|3.7|l that under Q* for 
each H in C 2 ([0,1]) 

M t (H) = Y t (H) - Y (H) - f ds Y S (AH) (3.9) 

Jo 

is a martingale with deterministic quadratic variation given by 

(M(H)) t = 2 f ds f du X {p{s,u)){VH) 2 . 
Jo Jo 

In particular, for each H, M t (H) is Brownian motion changed in time. 

Consider the semi-martingale Y s {T t - s H) for < s < t. Apply Ito's formula to 
derive equation l|3.5l) . with 



W t (H) - f dM s (T t - s H) 
Jo 



Wt{H) has a Gaussian distribution because the martingales M t (H) are Gaussian, 
being a deterministic time-change of a Brownian motion. The expression for the 
variance of Wt (H) follows from an elementary computation, as well as the fact that 
Wt{H) and Yq(G) are uncorrelated. This concludes the proof of the proposition. 

□ 

In view of (|3.5|l . to prove that Y N converges it remains to guarantee the conver- 
gence at the initial time: 

Proposition 3.3. Assume that Y N converges to a zero-mean Gaussian field Y 
with covariance denoted by <C ■, • 3>.' 

lim E.n[Y(H)Y{G)] = E[Y(H)Y(G)] =:<ff,G> . 

iV— »oo 



s 
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Then, Q converges to a generalized Ornstein- Uhlenbeck process with covariances 
given by 

E[Y t (H)Y s (G)} = <$:T t H,T s G-> 

r r 1 (3-io) 

+ 2/ dr X (p(r,u))(VT t _ r H)(u)(VT s - r G)(u) 
Jo Jo 

for allO <s <t< T, H, G m Cg([0, 1]). 

Proof. By Proposition 13.21 the sequence Q N is tight and all limit points satisfy 
(13.51) . Since W t and Y are zero- mean Gaussian random variables, so is Y t . To 
compute the covariance, it is enough to remind that the variables are uncorrelated. 
The first piece in formula l|3.10[) accounts for the covariance between Yo(T t H) and 
Y (T S G), while the last one for the covariance between W t {H) and W S (G). □ 

A noncquilibrium central limit theorem for the density field follows from the 
previous two results for processes starting from local Gibbs states. Indeed, fix a 
Lipschitz profile 7 : [0, 1] — > [0, 1] such that 7(0) = a, 7(1) = j3 and denote by 1/^ 

the product measure on {0, 1} Ajv associated to 7 so that 

v^Hx) = 1} = j(x/N) 

for x in A^r. In this case p N (x) = j(x/N), Lp N (x,y) — and p = 7. The first 
hypothesis in Ij3.2|l is satisfied because we assumed 7 to be Lipschitz. On the other 
hand, computing the characteristic functions of Y a N , it is easy to show (cf. 0) that 
Y N converges to a zero-mean Gaussian field with covariance given by 

E[Y{H)Y(G)] = [ X (l(u))H(u)G{u)du. 
Jo 

Therefore, by Propositions 13. 21 and 13.31 the density field converges to a generalized 
Ornstein-Uhlenbeck process: 

Corollary 3.4. Fix T > and a positive integer k > 5/2. Denote by Qn the 
probability measure on D([0, T], Ji-k) induced by the density fluctuation field Y N 
and the probability measure v ^i.y Then, Q N converges to the centered Gaussian 
probability measure Q with covariances given by 

E[Y t (H)Y s (G)} = f du X {l{u))T t H{u)T s G{u) 
Jo 

+ 2/ dr f du X (p(r,u))(VT t _ r H)(u)(VT s _ r G)(u) 
Jo Jo 

for all < s < t < T, H, G in Cg([0, 1]), where pt is the solution of the heat 
equation with initial condition 7. 

A similar result was obtained by De Masi et al. 0, for the one-dimensional 
symmetric exclusion process in infinite volume. This result was extended to higher 
dimensions by Ravishankar Chang and Yau [2] introduced a general method to 
prove non-stationary fluctuations of one-dimensional interacting particle systems. 

In Proposition 13.31 the asymptotic behavior of the covariance (|3.10ll as t f 00 
can be computed. Indeed, fix a function H in Cq([0, 1]) and set G = H, s = t. 
Since TtH vanishes as t ] 00, the first part of the covariance converges to 0. Since 
H vanishes at the boundary, since T t is the semigroup associated to the Laplacian 
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and since 2H\7H = Vif 2 , an integration by parts shows that the second part of 
the covariance (13.10(1 is equal to 

ds [ du X (p(s,u)){T t _ a H){u){d s T t _ s H){u) 
Jo 

t ,i 

ds / du[V X (p(s,uj)][V(T t _ s H(u)) 2 } . 
Jo 

Since 2Gd s G = d s G 2 , integrating by parts in time, since T t H vanishes in the limit 
1 1 oo and since the solution of the heat equation converges to the stationary profile 
p, the previous expression is equal to 

1 du X {p{u))H{uf - f ds [ du[d sX (p( S ,u))}(T t ^ s H(u)) 2 
Jo Jo 

t ,1 

ds / du[A X (p(s,u))}(T t _ s H(u)) 2 
Jo 

plus a term which vanishes in the limit. This sum is equal to 

du X (p{u))H{u) 2 - 2 f ds f du[Vp{s,u)] 2 {T t ^ s H{u) f 
Jo Jo 

because p is the solution of the heat equation. As t f oo, this expression converges 
to 

dux{p{u))H 2 {u) -{(3-a) 2 [ du H{u) {-h)- l H{u) 

Jo 

because Vp = (3 — a. We just recovered the covariance (|2.4|) of the density field 
under the stationary state. 

We turn now to the proof of Theorem 12. II Assume that the initial state p N is 
the stationary state g. We prove in Section 0] that the second condition in ((3.2(1 
is fulfilled. 

Fix k in K and recall the definition of the probability measure q introduced just 
before 12.4(1 . Let Q be the probability measure on C([0, T], H-k) corresponding to 
the stationary generalized Ornstein-Uhlenbeck process with mean and covariance 
given by 

E Q [Y t (H)Y s (G)] = E q [Y(T t _ s H)Y(G)] (3.11) 
for every < s < t and H, G in Ttk- 

Proposition 3.5. Fix T > and a positive integer k > 5/2. Denote by 
the probability measure on D([0, T],7i_fc) induced by the density fluctuation field 
Y N and the probability measure g . The sequence Qn converges weakly to the 
probability measure Q. 

Proof. Since v^r satisfies assumptions ((3.2(1 . ((3.3(1 . by Proposition ^. 21 Q N is tight 
and all limit points satisfy 1(3.5(1 . where Wt(H) is in this stationary context a zero- 
mean Gaussian variable with variance given by 

t ,i 

ds / du X (p(u))(T s VH(u)) 2 . 
Jo 

As t | oo, T t H vanishes in L 2 ([0, 1]). On the other hand, the computations per- 
formed just before the statement of this lemma show that the variance of Wt{H) 
converges to ((2.4(1 . Therefore, Wt(H) converges in distribution to a zero-mean 
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Gaussian variable with variance given by Ij2.4(l . Since the process is stationary, we 
just proved that the variables {Yt(H) : t > 0, H E Tik} have a zero mean Gaussian 
distribution with covariance given by (|2.4[1 . 

To compute the covariances E Q * [Y t (H)Y a (G)] = E Q * \Y t - a (H)Y (G)], < s < t 
it is enough to iterate relation (|3.9[1 to recover formula (|3.11|) . This concludes the 
proof of the lemma. □ 



We conclude this section proving that the sequence of probability measures Q n 
is tight and that all limit points are concentrated on continuous paths. 

To prove that the sequence Qn is tight we need to show that for every < t < T, 



lim limsupPpjv sup ||lt||-fc > A 



= 



and that 



lim lim sup P„j» 

S—0 N^oo 



w s (Y) > e 







for every e > 0. Here ws(Y) stands for the uniform modulus of continuity defined 

by 

w s (Y) = sup ||n-Y s ||_ fe 



\s-t\<6 
0<s,t<T 



We start with a key estimate. Recall the definition of the martingales M]' N (H), 
Mt' N {H) defined by ET7I) . 

Lemma 3.6. Fix a sequence of probability measures : N > 1} satisfying (I3.2[l . 
(I3.3() . There exists a finite constant C%, depending only on Co, such that for every 
.1 ■ 1- 



lim sup E, 



N- 



sup Y t {e 3 f < d/(l+T) 2 



0<t<T 



Proof. Recall and write Y t N (e 3 ) as M}' N (e 3 ) +Y N (e 3 ) + J * Y s N (A N ej)ds. We 
estimate these three terms separately. 

It follows from l|3.2[) that E^w [Yo(ej) 2 ] is bounded by a finite constant Ci, uni- 
formly in N and j. 



Since Ml' N (e 3 ) is a martingale, by Doob inequality, 



E„ 



sup \Mt- N (e 3 )f 

0<t<T 



< 41 



By definition of the martingale M^' N \e 3 ) and by l|3.9|) . the right hand side is equal 



to 



T 1 N-2 



ds-J2iVs(x + l)-Vs(x)] 2 {(V N e 3 )(x/N)Y + 0(A^ 



N 

By Lemma l3.1l as N f oo , this expression converges to 8 Jq dt x(/o(£, u))(Ve 3 >(u)) 2 du 
which is bounded by C\Tj 2 . 

Finally, by definition of e 3 and by Schwarz inequality, 



E„ 



sup ( / d S Y s (Ae 3 )) 2 } < Ci/TE^ 

0<t<T Wo 11 



dsY s {etf 



(3.12) 
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for some finite constant C\ . The previous expectation can be rewritten as 

(•<•)] 

i-T 

ds ■ 



f ds lt E ej{x/N) 2 p^(x)[l 

1 £ ejixMejiyMrffav) , 



where 

By Proposition ^. 41 



ft ( X ,V) = E u N [Vt(x);r] t (y)} 



sup sapN\(p?(x,y)\ < C x 

N>2 t>0 



(3.13) 



(3.14) 



for Ci = Co + (1/2)Cq. Hence, (|3.12|) is bounded above by CiT 2 j 4 , which concludes 
the proof of the lemma. □ 



Corollary 3.7. For k > 5/2, 

(a) limsupE^N 

N^oo 



sup ||lt||_ fc 
0<t<T 



< oo 



f&j lim limsupE^N sup y^Y t (ej) 2 j 2k 



0<t<T 







The proof of this result is similar to the one of Corollary XI. 3. 5 in [7] and therefore 
omitted. 

In view of Lemma f3. 61 and part (b) of Corollary 13. 71 in order to prove that the 
sequence Qn is tight, we only have to show that 



lim lim sup P„i\ 

s ~>0 N^oa 



sup {Y t {e j )-Y s {e j )} 2 >£ 

0<\s-t\<8 
0<s,t<T 



= 



for every j > 1 and e > 0. Fix j > 1 and recall the definition of the martingale 
Ml' N {e 3 ). Since Y t N {e 3 ) = Y N ( ej ) + M^ N {e ) + /* (e^ds, the previous state- 
ment follows from the next two claims: For every function G in Cq ([0, 1]) and every 
£ > 0, 



lim lim sup 

<5->0 at^oq 



lim lim sup I 



5^0 



N— >0O 



sup |M^(G)-M^(G)|>e 

0<s,t<T 
|t-«|<<5 

■i 



0. 



sup 

0<s,i<T 
|t-«l<<5 



dry r JV (A w G) 



> e 



The derivation of these estimates is similar to the proofs of Lemmata XI.3.7 and 
XI.3.8 in |2j if one keeps in mind the arguments presented in the proof of Lemma 
13.31 and the bound Ij3.14|l on the two point correlation function ip^ given by Ij3.13|l . 

4. Semidiscrete heat equation 

We prove in this section a bound on the two point correlation function tpf (x, y) 
introduced in (|3.13|l . Throughout this section, N > 2 is fixed. 

For the square of points C — {0, • ■ • , A^} 2 , consider the subsets V = {(x,y) £ 
C : < x < y < N} and its boundary dV = {(x, y)eC: x — or y — N}. Let 
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M = {/ : V U dV i— > R : f\ gv — 0} and denote by Ay the discrete Laplacian on 
M denned by 

(A$f)(x,y) = N 2 {f(x + l,y) + f(x-l,y) + f(x,y-l) + f(x,y + l)-4f(x,y)} 
if |x — y\ > 1 and 

(A^/)(x,x+l) = N 2 {f(x-l,x + l) + f(x,x + 2)-2f(x,x + l)} . 

Ay corresponds to the generator of a symmetric random walk on V U dV which is 
absorbed on dV. 

We start with an explicit formula for the total time spent by the random walk 
on the diagonal, which is expressed by the Green function or as the solution of the 
elliptic equation (— Ay (p N )(x, y) — CS y=x +i, C in K. Let (p N be the solution of 

{-A N ip N )(x, y) = C S y=x+ i for (x, y) in V , 

(4.1) 

tp N (x,y) = for (x,y) in dV . 
An elementary analysis shows that the unique solution of l|4.1|l is given by 

We turn now to maximum principles for solutions of homogeneous semidiscrete 
parabolic equations. Fix a function p N : A at — > M and let pf be the solution of 

d s pf (x) = A N p s (x) , x e A N 

p»(x)=p N (x), xeA N (4.3) 
p?(0) = a, p»(N) = /3, s>0. 

Lemma 4.1. Let p^ be the solution of (|4.3|) . Then, 

s>0 0<x<N — 1 0<:c<iv — 1 

Proof. Fix T > 0. Let 7t (a;) = p t (x + 1) - p t (a;) for < a: < JV — 1. Since 
dt1t(x) = (Ajv7t)(x) for 1 < x < A — 2, by the maximum principle, 

M = max sup |7 a (at)| = 

0<a;<W-l <s<T 

maxi max \p N (x + 1) - p N (x)\ , sup |p t (l)-a|, sup |/3 - p t (A - 1)| I . 

0<a;<W-l 0<t<T 0<t<T J 

We claim that the maximum is attained at t = 0. To show this assume, without 
loss of generality, that there exists to £ (0, T] such that M = |pt (l) — ck|. By (|4.3|) 
with i = lwe have, 

fl a p.(l) = -N 2 p s (l) + N 2 a + N 2 (p s (2) - p,(l)) 

for any < s < T. Thus, multiplying by e sN , grouping the terms conveniently 
and integrating on [Q,i] we get that 

p t (l) = e- tN2 p (l) + a(l - e~ tN2 ) + f e^'^ N 2 (p s (2) - p s (l))ds 



so that 

Pt(l) - a = e- tN2 (p (l) - a) + I A 2 (p s (2) - Ps (l))ds 
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Using the assumption made on |pt (l) — a\, we deduce from this identity that 



M = K(l)-a| < e- toAr -| P0 (l)-a| + (1 - )M , 

which reduces to M < |po(l) — ct\. This concludes the proof of the lemma. □ 

Fix a function h : V — > K and denote by f s the solution of the semidiscrete heat 
equation 

d s f s (x,y) = A$f s (x,y) ,(x,y) eV , 

fo{x,y) = h(x,y) , (x,y) E V (4.4) 
f s (x,y) = 0, 0<s<T, (x,y)edV. 

Lemma 4.2. Let f satisfy j4-4\ l- Then, the maximum value of f on [0, T] x (VUdV) 
is attained at a point (to,xo,yo) such that to — or (xq, yo) £ dV . 

Proof. The proof is the same as that of the maximum principle for the usual heat 
equation. It uses that if the maximum is attained at an interior point (to,xo,yo) 
then d to ft {xo,yo) > and A^f to {x Q ,y ) < 0. □ 

Fix a function h : V — > K and a function g : M + xK^l. Consider the following 
nonhomogeneous parabolic equation, 

d s fs(x, y) = A$ip s (x, y) + g s (x, y) , (x, y) € V , 
ip (x, y) = h(x, y) , (x,y) G V , (4.5) 
<p s (x,y) = 0, 0< S <T, {x,y)edV, 

Denote by || • ||ioo(y) the sup norm: ||/i||;oo(y) = max ze y 

Lemma 4.3. Fix T > and assume that the function g s is supported on the line 
y = x + 1 . Then, 

sup H^ll^y) < ||/i||i°o (v) + . 1 . sup max \g t {x,x + l)\. 

0<t<T 4 W V — tj O^KT 1 ^- 2 

Proof. It is not difficult to see that Ay is a symmetric, negative operator on M. = 
{/ : V U dV i — ► M : f\ dv = 0}. It generates in particular a strongly continuous 

semigroup {e sA ™ : s > 0} on M. and the solution <p s of (|4.5|) can be written in the 
form 

<p s = e sA $h + f eS s -^g t dt. 



(i 



Since e sAv h is the solution of l|4.4(l . by the maximum principle stated in Lemma 

E21 

||e sA "^||^ (y) < \\h\\ lx(v) 
and e v (x, y) > 0. On the other hand, since gt is supported on the line y = x + 1, 

/ {e^^JOcy)^ < / Ve( s -*' A "(x,y,x',x' + l)| 5t (x',^ + l)|^ 
./o x , =1 

< / e tA ™ (x,y,x' ,x' + l)dt sup max \gt(x', x' + 1)| 
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provided s < T. To conclude the proof of the lemma it remains to show that 
the integral is bounded by (4[iV — 1]) _ , uniformly in x, y and s. This integral is 
bounded above by 

v N (x,y) = / J2e tA v(x,y,x',x' + l)dt (4.6) 

Jo x'=l 

which satisfies the equation Ay v^(x,y) = —6y= x +i. By l|4.2J) . 

1 x 



v N (x,y) = — 



and this expression is less than or equal to (4[JV — 1]) 1 because x < y. This 
concludes the proof of the lemma. □ 



We are now in a position to state the main result of this section. 

Proposition 4.4. Consider a probability measure fi N on {0, 1} Ajv satisfying (|3.2[1 . 
Denote by tp^ the two point correlation function: 



<p?(x,y) = {r) t (x)-p»(x)}{r, t (y)-p?(y)} 



wh 



N 

ere p t is 



the solution of (|4.3|) . Then, 



supii^ri 

t>0 



< 



2C Q 



2N 



for all N>2. 



Proof. A simple computation shows that ip^(x,y) is the solution of (|4.5|l with 
h = £yv [77(3;); 77(2/)] and g t {x,y) = -(V ' NP?){x) 2 8 y=x+ i. By Lemma EH] and 
assumption (|3.2|) . g t , which is supported on the line y — x + 1, is absolutely bounded 
by Cq. Therefore, by Lemma 14.31 iVsup t > Hyj^ ||;oo(y) is less than or equal to 
Co + Cq/2 since N > 2. This concludes the proof of the lemma. □ 



It remains to check that stationary state v£p satisfies the assumptions of the 
previous proposition. Recall the definition of p N (x) given just before l|2.2|l and 
denote by (p N (x, y), l<x<y<N — 1 the two point correlation function: 



V "(x,y) = E <p [Mx)-p N (x)}{T,{y)-p N (y)} 



(4.7) 



Since E v n [L n {t)(x) - p N {x)} {q(y) - p N (y)}} = for all 1 < x < y < N - 1, we 

obtain that (p N (x,y) is the solution of the discrete differential equation l|4.1() with 
C — (3 — a. Therefore, 



(/3 -a) 2 x 



N 



'■T_ ( 1 _ y_\ 

1 TV V NJ 



In particular, (|3.2|) is satisfied and we may apply Proposition 14.41 
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